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ABSTRACT 

High-resolution infrared spectroscopy in the 2.3-4.6 /im region is reported 
for the pecuhar A supergiant, single-hned spectroscopic binary HR 4049. Lines 
from the CO fundamental and first overtone, OH fundamental, and several H2O 
vibration-rotation transitions have been observed in the near-infrared spectrum. 
The spectrum of HR 4049 appears principally in emission through the 3 and 4.6 
/im region and in absorption in the 2 //m region. The 4.6 /Um spectrum shows a 
rich 'forest' of emission lines. All the spectral lines observed in the 2.3-4.6 /im 
spectrum are shown to be circumbinary in origin. The presence of OH and H2O 
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lines confirm the oxygen-rich nature of the circumbinary gas which is in contrast 
to the previously detected carbon-rich materiaL The emission and absorption 
fine profiles show that the circumbinary gas is located in a thin, rotating layer 
near the dust disk. The properties of the dust and gas circumbinary disk and 
the spectroscopic orbit yield masses for the individual stars, IsAa i ~0.58 M0 and 
Mm V ~0.34 Mq. Gas in the disk also has an outward fiow with a velocity of ^1 
km s~^. The severe depletion of refractory elements but near-solar abundances 
of volatile elements observed in HR 4049 results from abundance winnowing. 
The separation of the volatiles from the grains in the disk and the subsequent 
accretion by the star are discussed. Contrary to prior reports, the HR 4049 
carbon and oxygen isotopic abundances are typical AGB values: ^^C/^^C=6^4 
and 160/17O200. 

Subject headings: accretion disks — stars:abundances — stars:AGB and post- 
AGB — stars: chemically peculiar — stars:evolution — stars:winds,outfiows 



Introduction 



HR 4049 is the pr ototype for a class of pe culiar, post- AGB, single-lined, long-period, 
spectroscopic binaries (Ivan Winckel et al.lll995l ). The primary star of these binaries is an 
early-type supergiant with very peculiar abundances. Most objects in this class exhibit 
strong infrared excesses of circumstellar origin with carbon-rich circumstellar dust typically 
present. The combination of carbon-rich material, high luminosity, and location out of 
the galactic plane forms the basis for the post-AGB designation. The peculiar designation 
stems from a photospheric abundance pattern characterized by a severe deficiency of refrac- 
tory (high dust condensation temperature) elements and a near-solar abundance for volatile 
(low dust condensation temperature) elements. The abundance anomalies indicate that the 
present photosphere contains material from which refractory elements have been very largely 
removed, i.e., a winnowing of dust from gas has occurred. 

The basic characteristics of the prototype object HR 4049 are well known. The photo- 
sphere has an effective temperature of about 7500 K and shows extreme abundance differ- 
ences betwe en refractory and volat i le elements: for exa mple, HR 4049 has [Fe/H]~-4.8 but 
[S/H]~-0. 2 JWaelkens et"aDll991bl : lTakada-Hidailll99oh . The orbital period of HR 4049 is 
430 dav s (IWaelkens et al.lll991al) . leading to a minimum separation between the two stars of 



190 R0. iBakker et al. 



Jl998h pointed out that the orbit requires a phase of common envelope 
evolution when the primary was at the tip of the AGB and had a radius of ~ 250 R©. This 
phase altered the masses and abundances of the components. 
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The HR 4049 infrared excess is pronounced redward of ~1 fim and is very well fit by a 
single blackbody at a temperature of about 1150 K a nd attributed to radi ation from the tall 
inner walls of an optically thick circumbinary disk (jPominik et al.ll2003l ). A circumbinary 
Keplerian rotating disk appears a common feature of the HR 4049 class of post-AGB binaries 
( jPe Ruyter et al.l 120061 ) . The inner walls of the HR 4049 circumbinary disk are ~ 10 AU 
from the binary or 50 times the radius of the supergiant. Our line of sight to the binary 
nearly grazes the edge of the disk; the angle of inclination of the line of sight to the normal 
to the disk is about 60°. A cartoon of the system is shown in Figure 1. 



Superimposed on the infrared dust continuum are emission features. IWaters et al. 



( 119891 ) detected features due to polycycli c aromatic hydr o carbo n s (PAHs). Thi s resul t was 
confirmed by the ISO/SWS spectrum of iBeintema et al.l (119961 ). iGeballe et al.l (119891 ) con- 
firmed a C-rich circumstellar environment by detecting 3.43 and 3.53 / im emission feature s 
later identifi ed with hydrogen-terrn inated crystalline facets of diamond (IGuillois et al.lll999l ). 
Remarkably, iDominik et al.l (120031 ) report that the gas molecular species seen in the infrared 
(ISO) spectrum are those expected of an 0-rich mixture, suggesting that the HR 4049 cir- 
cumbinary environment is a blend of C-rich dust and 0-rich gas. 



Optical spectroscopy offers some information on the disk's gas. Bakker et al.l (19981) re- 
port c hanges in the Ha line profile with orbital phase. iBakker et al.l (119981 ) and lBakker et al. 
( 119961 ) detected a broad (~15 km s~^) stationary emission component in the Na I D2 and [O 



I] 6300 A lines which they attributed to the circumbinary disk. However, the infrared offers 



much more readily interpreted signatures of gas in and around the binary. iLambert et al. 
( 119881 ) detected the ^^C^ ^O first overtone sp e ctrum in absorption with an excitation tem 



perature of about 300 K. ICami fc Yamamural (I2OOII ) on analyzing an ISO/SWS spectrum of 
CO2 emission bands found strong contributions from isotopomers containing ^''O and ^^O 
which they interpreted as ^^O/^^O = 8.3 ± 2.3 and ^^O/^^O = 6.9±0.9. 

The origins of the HR 4049 class of chemically peculiar supergiants and the structure of 
their circumbinary /circumstellar material remain ill-understood. New observational attacks 
appear to be essential. In this paper we report on a detailed look at several regions of the 
2-5 ^m infrared spectrum of HR 4049. 



2. Observations & Data Reduction 



The spectrum HR 4049 was observed at high resolution at a number of near-infrared 
wavelengths in the 2.3 to 4.6 /xm region using the 8 m G emini South telescope and th e 
NOAO high- resolution near- infrared Phoenix spectrometer ( iHinkle et al.lll998l . I2OOOI . l2003l ). 
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Phoenix is a cryogenically cooled echelle spectrograph that uses order separating fihers to 
isolate sections of individual echelle orders. The detector is a 1024x1024 InSb Aladdin II 
array. Phoenix is not cross dispersed and the size of the detector in the dispersion direction 
limits the wavelength coverage in a single exposure to about 0.5%, i.e. 1550 km s~^, which 
is 0.012 /im at 2.3 (22 cm"^ at 4300 cm'^) and 0.024 /xm at 4.6 ^m (11 cm"^ at 2100 
cm^^). One edge of the detector is blemished so the wavelength coverage is typically trimmed 
a few percent to avoid this area. Wavelength coverage is limited overall to 0.9-5.5 /im by 
the InSb detector material. All the spectra discussed here were observed with the widest 
(0.35 arcsecond) slit resulting in a spectral resolution of R=A/AA= 50,000. The central 
wavelengths of the regions observed are listed in Table 1. 

The thermal brilliance of the sky makes observations longward of ~4 fim much more 
difficult than in the non-thermal 1-2.4/xm region. However, for a bright star like HR 4049 this 
only slightly increases the already short integration t ime. Therrn al infrared observations were 



done using standard infrared observing techniques (jJoycdll992l ). Each observation consists 
of multiple integrations at several different positions along the slit, typically separated by 
4" on the sky. At thermal infrared wavelengths the telluric lines are in emission. In order not 
to saturate the telluric emission lines the limiting exposure time is about 30 seconds at 4.6 
//m. For longer integration times, multiple exposures can be coadded in the array controller 
to make up a single exposure. However, HR 4049 is so bright that total exposure times of 
only 10 to 20 seconds were required. The delivered image FWHM at the spectrograph varied 
from 0.25"- 0.80"during the nights that spectra were taken. With the positions along the 
slit separated by several arcseconds the resulting spectral images were well separated on the 
detector. 

An average flat observation minus an average dark observation was divided into each 
frame observed and frames with the star at differentplaces along the slit were then differenced 
and the spectrum extracted using standard IRAF]^ routines. A hot star, with no intrinsic 
spectral lines in the regions observed, was also observed at each wavelength setting. The 
hot star was observed at airmass near that of HR 4049 and the HR 4049 spectrum was 
later divided by the hot star spectrum to ratio the telluric spectrum from the HR 4049 star 
spectrum. Wavelength calibrations were computed by using a set of telluric wavelengths 
obtained from the hot star spectra. The wavelength calibration yielded residuals of typically 
0.25 km s~\ 

Observations of HR 4049 were taken in the 2 and 3 yum region as well as the 4.6 /im 



^The IRAF software is distributed by the National Optical Astronomy Observatories under contract with 
the National Science Foundation. 
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region. For a bright star the stellar signal is much stronger than background radiation in 
these spectral regions and as a result the observations are much less challenging than 4.6 fim 
observations . The Phoeii i x obs erving technique for this spectral region has previously been 
described in ISmith et al.l (120021 ) . 



3. Analysis of the Spectra 

Observations of the 2.3, 3.0, and 4.6 /im regions reveal lines from just three molecular 
species: CO, OH, and H2O. As discussed in §1, HR 4049 has a carbon-rich circumstellar 
envelope but we did not identify any molecules associated with conditions where C>0. In 
the 4.6 ^m region we searched for C3 and CN which should be prominent if C>0. Rather, 
we report on the new detection of a rich 4.6 /um forest of CO fundamental and H2O lines. 
The 4.6 fim region atomic hydrogen lines Pfund /3 and Humphreys e, if present, are blended 
with H2O lines. The 4.6 fim HR 4049 emission line spectrum is very rich making the identi- 
fication of occasional atomic features problematic. The circumstellar continu um at 4.6 iim is 



appro ximately 25 times more intense than the continuum of the supergiant (jPominik et al. 



20031 ) so features of stellar origin will be highly veiled. At 3 yum we make a first detection of 
the OH fundamental vibration- rotation lines in HR 4049. Before exploring the emission line 
spectrum, we revisit in §3.1 the 2.3 /im CO first overtone spectrum detected previously by 



Lambert et al.l (119881 ). 



All velocities in this paper are heliocentric. In order to compare heliocentric velocities 
of HR 4049 with microwave observations add -11.6 km s~^ to convert to the local standard 
of rest. 



3.1. CO First Overtone 



Our observations confirm a nd extend the discove ry of first-overtone (Av=2) vibration- 
rotation CO lines in absorption (ILambert et al.lll988l ). At the 2.3 /im wavelength of the CO 
first overtone, the continuum from the dust is about four times that of the supergiant. Thus, 
the CO absorption lines should be formed along the lines of sight to the circumbinary disk. 
By inspection (Figure 2), it is apparent that the rotational and vibrational temperatures 
are low; high rotational lines of the 2-0 band are absent and the 3-1 (and higher) bands are 
weak or absent relative to the low rotational lines of the 2-0 band. All of the prominent lines 
are attributable to the most common isotopic variant, ^^C^^O, but weak 2-0 ^^CO lines are 
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detectabl^. 

The first set of observations from February 2002 show the CO lines at a radial velocity 
of -33±0.5 km s~^. The lines reach maximum strength at J"~7, suggesting a rotational 
excitation temperature of ~300 K. At J" =7, the R branch lines are about 17% deep and are 
resolved with a FWHM of 16 km s~^ compared to the instrumental resolution of 6 km s~^. 
Weak ^^CO lines are detected with depths for the R18 to R23 lines of about 4%. Comparing 
lines of similar excitation suggests that ^^C/^^C~10. The observed regions cover the strongest 
predicted C^^O lines (2-0 lines near J"~ 7) but these lines can not be convincingly identified 
in the spectra (Figure 2) demanding C^^O/C^^O ^100. A yet more stringent limit can be 
applied (§4) by modeling. 

After the original observations, additional data were collected to extend the excitation 
range of the lines. For instance, observations made in December 2002 included higher J 
lines than observed previously. Observations in December 2005 covered the low J P branch 
required for curve-of-growth analysis. Some wavelength intervals were reobserved over the 
2002 - 2005 interval to check for variability. The radial velocity was in all cases unchanged 
from the -33 km s~^ measured in the original data set. This veloci ty is nearly equal t o the 



-32.09 ±0.13 km s ^ systemic velocity of the spectroscopic binary (IBakker et al.lll998l ). 



The CO first overtone line profiles are symmetrical with no hint of an emission compo- 
nent. The line strengths sho wed no tempora l varia bility. In fact, the line intensities are quite 



similar to those reported by lLambert et al.l (119881 ). Similarly, the velocity is identical to the 



earlier reported value. The profile of the lowest excitation line, 2-0 RO, differs from others in 
that it possibly has a weak blue-shifted component. However, the 2-0 RO line lies in a region 
with a complex telluric spectrum which is difficult to ratio out of the HR 4049 spectrum. In 
the February 2003 spectrum, the 2-0 RO line appears to have components at -23.6 and -32.9 
km s~^. On other dates, the blue-shifted component is less clearly resolved suggesting that 
it is either of variable velocity, affected by overlying emission of variable intensity and/or 
velocity, or a relic of the reduction process. If this blue-shifted component exists, it must 
originate in very cold gas (T ^ 5 K) because the component is not detectable in the 2-0 Rl 
line. 



^We follow the convention of omitting the superscript mass number for the most common isotope. Hence 
i2(-;16q appears as CO, etc. 
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3.2. CO Fundamental 

In sharp contrast to the spectrum at 2.3 where a sparse collection of weak absorption 
lines are found, the spectrum near 4.6 fim is rich in emission lines (Figure 3). Emission 
lines were identified from four isotopic variants: CO, ^^CO, C^''0, and C^^O with roughly 
equal intensities for all variants. The rarer isotopic variants ^^C^'^'O, ^^C^^O, and ^^CO were 
searched for but are not present. Absorption below the local continuum is also seen in the 
profiles of the lowest excitation 1-0 CO lines in this interval. The observed spectral interval 
provides lines mostly from the 1-0 and 2-1 CO bands but a few R branch lines of the 3-2 
CO band are clearly present. The maximum observable rotational level, J"~30, is similar to 
that seen in the CO first overtone. Table 2 lists the detected fundamental lines of the four 
CO isotopic forms. 

Blending with other lines of different CO isotopes, vibration-rotation transitions, or H2O 
lines is common and results in an apparent variety of profiles. All unblended emission lines 
arc double-peaked - see Figures 4-8. For all but the lowest excitation lines the blue and 
red peaks are of similar intensity, but characteristically with the blue peak slightly weaker 
than the red, and occur at velocities of -38.9 ±0.4 and -28.4 ± 0.5 km s~\ respectively. The 
central valley of the emission profile has a velocity of -33.7 ± 0.3 km s~^. Velocities are not 
dependent on the isotopic species. The observed, i.e. uncorrected for instrumental profile, 
full- width at zero intensity (FWZI) of the emission profile for the weaker lines is ~ 27 km 
s~^, with stronger lines having FWZI up to ~ 35 km s~^. Due to the high line density, 
the FWZI is a difficult parameter to measure and our values carry an uncertainty of several 
km s~^. The observed full-width at half maximum (FWHM) similarly depends on the line 
strength but much less dramatically than the FWZI. Typical FWHM values are ~ 19 km 

At the observed resolution of A/AA=50000, the instrument profile has a significant 
impact on the observed line profiles and FWZI. Assuming a Gaussian instrumental profile 
equal to the 6 km s~^ FWHM spectral resolution, deconvolution of this instrumental profile 
from the observed line profile gives a true FWZI of 18 km s~^. The line profiles are strongly 
smoothed by the instrumental profile. The observations show blue and red sides of the CO 
fines rising ~20% above the body of typical CO emission lines with the peak at each edge 
having a FWHM of ~ 4 km s~^. For more strongly saturated fines, e.g. low excitation ^^CO 
fines the FWHM of the blue and red emission spikes are ~8 km s~^. The intrinsic profile 
clearly has much stronger emission peaks. 

The combined absorption-emission profile for the 1-0 CO lines is shown best by the R2 
and R5 fines (Figure 9). After allowance for similar blends, the profiles of the CO Rl, 2, 3, 
and 4 fines can be judged very similar to that of the R5 line. The profile of the RO line is 
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possibly of the same type but blending is more severe. The R5 profile is almost a P Cygni 
profile: blue absorption accompanied by red emission. However, the absorption component 
is not strongly blue shifted but has a velocity very similar to that of the absorption lines in 
the other fundamental lines and to that of the 2-0 lines. 

Absorption below the continuum is seen only in the 1-0 P and R branch CO lines (Figure 
3 - 8). The observed interval includes the PI, P2, P3 and RO, Rl, R2, R3, R4, and R5 lines 
with definite absorption below the local continuum seen in all these lines except PI and RO. 
The strength of the absorption at R5 suggests absorption below the local continuum should 
be detectable to higher J hues of the R branch. Due to the isotopic shifts, the low J 1-0 
lines of the isotopic variants are not in the observed interval. The lowest member of the R 
branch in our spectra is J" = 9 for ^^CO, J" = 10 for C^^O, and J" = 3 for C^^O. All the 
1-0 lines regardless of isotopic species have a stronger central valley than the vibrationally 
excited transitions. The valley almost reaches the local continuum for the 1-0 ^^CO lines 
(note i^CO 1-0 RIO and Rll in Figure 7). 

The central valley in the line profiles becomes asymmetric for the stronger lines. Self- 
absorption of the blue emission is obvious when comparing the strength of the blue and red 
emission in profiles of the 1-0 CO lines (Figure 9). The 1-0 CO line central absorptions 
are on average 76% broader on the blue side than the red side. The weaker ^^CO 1-0 
central absorptions are 20% broader on the blue side . The central absor ption is systemically 



blue shifted relative to the 7 velocity of the binary (iBakker et al.lll998l ) for the very lowest 
excitation lines. The shift increases with decreasing J", with a shift of 1 km s~^ for R5 and 
3.5 km s~^ for R2 (Figure 10). 

The intensities of the emission lines of CO, ^^CO, C^^O, and C^'^0 are remarkably similar 
and quite different from the abundance ratios estimated from the 2-0 lines. Peak intensities 
of the following representative unblended lines illustrate this point: 
CO: 1-0 R2 28%, 2-1 Rll 27%, 3-2 Rll 11% 
i^CO: 1-0 R12 24%, 2-1 R17 14% 
C^^O: 1-0 Rll 9%, 2-1 R17 14% 
C^^O: 1-0 R14 14% 



In contrast to the ratio CO/^^CO ~ 10 from the 2-0 lines, the CO/^^CO intensity ratio 
from 1-0 and 2-1 lines of similar J is about 1.5. Even more striking is the appearance of 
fundamental lines of C^^O and C^^O with intensities about one-half that of similar lines of 
CO. Yet, CO/C^^O > 100 from the first-overtone lines. The simplest interpretation of these 
contrasting ratios is that emitting regions are optically thick in all the observed fundamental 
lines. As is well known, the first-overtone transitions are much weaker than the fundamental 



- 9 - 



lines. Optically thin emission in fundamental and first-overtone lines from a common upper 
state in the second vibrational level will differ by a factor of about 100 in flux. In the case of 
absorption from a common state in the ground vibrational level, the absorption coefficient 
of the 1-0 line is similarly about factor of 100 stronger than the 2-0 line. 

The strengthening of the central absorption for the lowest energy vibrational transition, 
the asymmetric absorption, and absorption below the continuum require the presence of an 
absorbing gas cooler than the emitting gas. This absorbing gas has a velocity shifted to the 
blue of the system barycentric velocity (-32.1 km s~^) by 1 to 3.5 km s~^. The emitting gas 
covers a ~18 km s~^ range of velocity but is also shifted by ~-1.5 km s~^ relative to the 
barycentric velocity. 



3.3. OH Fundamental 

The lowest excitation OH vibration-rotation 1-0 lines are in a region of considerable 
telluric obscuration. J" =4.5 is the lowest OH level accessible under typical water vapor 
conditions (a few mm of precipitable H2O) at Gemini South. However, a suitable order 
sorting filter was not available for the J" =4.5 wavelength. An observation was made of the 
P branch line region for J" =5.5. The ^H OH ground state results in A-doubled rotational 
levels, so each rotational line is divided into four components. As a result, in spite of the 
large rotational line spacing for OH, several OH lines can appear in a Phoenix spectrum 
taken with a single grating setting. The 3.0 /xm 1-0 P2/5.5 and P2e5.5 lines were detected in 
the spectrum of HR 4049 (Fig. 11). This spectral region has considerable telluric absorption. 
The removal of this absorption results in variable noise in the ratioed spectrum. 

The OH lines the CO lines, seen in emission. The profiles are similar to those 

of CO with double peaked profiles of observed FWZI ~26 km s~^ centered at -35 km s~^. 
The two OH lines observed are just 5% above the continuum. These were the only lines that 
were detected in the 3.0 /im spectral region observed. 



3.4. H2O Vibration-Rotation Lines 

The asymmetric top molecule H2O is known for the complexity, apparent lack of rota- 
tional structure, and richness of its spectrum. As a result H2O lines are mu ch more challeng- 



ing t o identify than vibration-rotation lines of simple diatomic molecules (IHinkle &: Barnes 



19791 ). Emission lines are clearly present in the 4.6 fim HR 4049 spectrum from three 



vibration-rotation bands: U2, fi — 1^2, and ^3 — ^2- With the above caveats on the H2O 
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spectrum and based on the tentative identification of four lines, tlie vibrationally excited 
band 2^2 — possibly also contributes to the 4.6 /im spectrum. The strongest observed 
H2O transitions are from the — band. A number of lines identified with this band have 
intensities ~20% above continuum. Typical H2O lines are weaker than typical CO lines, 
with many of the H2O lines identified having intensity <10% above continuum. 

Table 3 presents a list of the H2O lines tentatively identified in the HR 4049 spectrum. 
In Figures 4-9 these lines are labeled on the spectrum of HR 4049. Many lines (e.g. Figure 
9) are unblended and clearly present. However, a fairly large number are blended with CO 
or other H2O lines. Due to the overlapping H2O energy levels, the line strengths of H2O 
lines can vary significantly between adjacent vibration-rotational transitions and, hence, the 
contribution of a H2O line to a blend is uncertain. 

The band strength, 5*°, is a factor of five lower for the vi — band than the 1^3 — V2 
band. However, the band strength for th e V2 band is more than 10^ higher than that of either 



of these bands (IMcClatchey et al.lll973l ). The origin of the V2 band is ~1.5 [xm red of the 
region observed. While it would be of interest to observe the lowest excitation V2 lines, for 
ground based observers the telluric V2 lines are very strong and prevent observations in 6 /im 
region. The 2z/2 — V2 band has similar band strength to the v\ — V2 and 2/3 — V2 combination 
bands but, unlike these bands which have origins in the regions observed, 2v2 — V2 has an 
origin near that of V2- This adds to our suspicion of the 2z/2 — V2 identifications. 

Like CO and OH lines the 4.6 /im H2O lines have a double peaked profile with emission 
peaks at -38.6 and -29.2 km s~^ and absorption at -33.6 km s~^. The observed FWZI of the 
weaker H2O lines is ~25 km s~^, perhaps slightly more narrow than the CO lines. 

In addition to the observed 4.6 /xm transitions, H2O also has low excitation transitions 
in the 3.0 and 2.3 /xm regions. In particular the 1^3 band crosses the 3 /im region and has a 
band strength similar to that of the 1^2 . The v\ band is also present in the 3 /xm region and, 
while weaker than V2 or z/3, is a much stronger transition than the combination bands seen 
at 4.6 /im. Our 3.0 /im observation has an uneven continuum perhaps as a result of weak 
emission features. We undertook a detailed search for H2O lines but failed to identify any 
3.0 /im H2O lines. Future searches of this region for H2O lines using higher signal-to-noise 
data and wider wavelength coverage are justified. On the other hand, our spectra in the 2.3 
/im region are of very high quality with broad wavelength coverage and this region is clearly 
devoid of any contribution from H2O. 



3.5. Line Profile Overview 



In summary of the above subsections, CO vibration-rotation fundamental lines of four 
differently isotopically substituted species and H2O vibration-rotation lines populate the 4.6 
yum region. These lines all have double peaked emission profiles of total (including both 
peaks) FWHM ~ 19 km s~^. There is little difference of intensity between lines of different 
isotopes. The lowest excitation fines, which are only seen in ^^C^^O in the wavelength range 
observed, have a central absorption as much as 20% below the local continuum. This central 
absorption overwhelms the bhiest of the double peaks in the emission profile but the extreme 
bluest edge of the emission remains. Examples of observed CO fundamental and H2O line 
profiles are given in Fig. 9. OH lines from the fundamental vibration-rotation transition were 
seen in the 3 /xm region. The OH lines are in emission with double peaked profiles similar to 
those seen in the CO fundamental and H2O lines. The CO vibration- rotation first-overtone 
transition appears in absorption in the 2.3 fxm region. The absorption lines are nearly as 
broad as the CO emission lines, FWHM ~ 16 km s"^, but have simple Gaussian profiles and 
exhibit a range of fine depths suggesting the fines are optically thin. ^^CO dominates but 
weaker -"^^CO is detectable. The oxygen isotopes are not present. All spectral lines in the 2-5 
fj,m region have a small 1 km s~^) shift blue of the systemic velocity. 



4. Modeling the Molecular Probes 

The fundamental spectrum presents a difficult analysis task. The CO is seen in both 
emission and, for the very lowest excitation CO lines, absorption. The small change in 
intensity for emission lines over the full range of isotopes and over a large range of excitation 
energy (rotational levels J"=0 to 30 and vibrational transitions 1-0 to 3-2) clearly indicates 
that the emission lines are very saturated. A detailed investigation of these strongly saturated 
lines would require detailed radiative transfer and disk modeling beyond the scope of this 
paper. 

However, analysis of the CO first overtone lines, which are seen in absorption, is a much 
more tractable problem. The observations cover nearly all ^^C^^O 2-0 R branch lines from 
J"=0 through the highest detectable R branch line at J" =35. The largest interval of the 2-0 
R branch not observed is R24 through R28. The 2-0 P branch was observed from PI through 
P8. The 3-1 R branch from J" > 4 also fies in the observed region. Equivalent widths of the 
first overtone CO line profiles were measured from the fully processed, normalized spectra. 
The lines were measured both by summing the absorption area and second by Gaussian 
fits to the line profiles. Uncertainties were estimated from the mean deviations from the 
Gaussian fits and by the formal uncertainty in the fitted continuum level. 
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Equivalent width data was used to produce an excitation plot (Figure 12). The log- 
linear increase of line strength with excitation level demonstrates that the high-J v=2-0 lines 
(J>15) are optically thin (i.e. r < 0.7). A least squares fit to the excitation plot of these 
data requires a temperature of ~550 K. However, the fit to these lines underestimates the 
column density of the low-J lines. To infer the temperature and optical depth of the low-J 
lines, we extrapolated the column density of the hot gas (inferred from the high-J lines) and 
subtracted that from the measured column density. 

In order to correct for the effects of saturation, column de nsities and level populations 



were determined from a curve of growth (COG) analysis (c.f. ISpitzerl Il978l : iBrittain et al. 



20051 ) . which relates the measured equivalent widths to column densities by taking into 
account the effects of opacity on a Gaussian line profile. The derived column density for a 
measured equivalent width only depends upon one parameter, the Doppler broadening of the 
line, b = cr/jM5/l-665, where aRMs is the RMS linewidth. To find the value of b, we apply two 
complementary methods: comparison of the P and R branch lines and the linearization of 
the excitation plot. A key assumption is that the small scale line broadening results entirely 
from thermal broadening. The resolved line profiles which are seen in the spectra indicate 
an additional large scale broadening mechanism which will be discussed in §5.1. 

CO exhibits absorption lines in both P (J"=J'+1) and R (J"=J'-1) branches, which 
have different oscillator strengths yet probe the same energy levels, e.g., the PI absorption 
line originates from the same J=l level as the Rl absorption line. Any differences in the 
column density derived from lines that share a common level must be due to optical depth, 
which is related to b. The line width, b, can be used to determine the optical depth and 
adjusted so that the derived level populations from the two branches agree as closely as 
possible. 

The (v=0, low-J) transitions are thermalized at densities as low as ylh ~ 10'^"'' cm~^, 
and at even lower densities due to radiative trapping in the rotational lines with high opacity. 
Therefore, the low-J lines are the ones most likely to exhibit a thermal population distri- 
bution. The line width that best linearizes a plot of the level populations to a common 
temperature in an excitation diagram is used. 

Subject to the above constraints, the best fitting b value in the COG analysis for HR 
4049 is 0.5±0.1 km s^^. The consistency of all data to this common velocity dispersion is 
depicted in the excitation plot of Figure 12. With a measurement of b, equivalent widths 
can be directly related to column density. The column density from fitting the 2-0 'high-J' 
lines is 4.6±0.3x 10^'' cm~^ at a temperature of 530±20 K. The column density of the 'low-J' 
lines is 1.6±0.2xl0^^ cm~^ at a temperature of 40±10 K. Uncertainty in the hot N(^^CO) 
from the overtone lines, estimated from the measurement of unsaturated lines, is small and 
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dominated by measurement errors in the equivalent widths of the hues. The uncertainty in 
the cold gas is dominated by the uncertainty in b. Assuming that the 0.5 km s~^ b value 
for the cold gas apphes to all the spectral lines, the opacity of the most optically thick line 
is ~1.5. Increasing the b value for the hot gas lowers the optical thickness of the higher 
excitation lines. 

The detection of weak 3-1 lines allows a check on vibrational LTE in the gas. The 
column of CO in the v=l state (from the v=3-l lines) is (5.7±1.3)xl0^^ and the temperature 
is 540±80K. This is consistent with the temperature for the hot component of the v=0 ^^CO 
and ^'^CO branches (530±20 and 570±40 respectively). The combined 2-0 and 3-1 data 
give a rotational temperature of 620±20K. The vibrational temperature is 700±50K. This 
is consistent with a slight overpopulation of the v=l state although the relative rotational 
populations are consistent. The vibrational temperature is more uncertain than the other 
temperatures and evidence for non-LTE populations is weak. 

Using the best fitting b value, the column density can be determined for other isotopic 
lines in the spectrum. The corresponding column density of ^^CO is 2.3±0.3xl0^^ cm~^ at a 
temperature of 570±40 K. Comparing lines of similar excitation, ^^C/^^C ratio is First 
overtone C^^O lines could not be detected. The strongest lines, assuming a 550 K excitation 
temperature, that are clear of both major telluric features and blending CO lines are R5 
and R8 (Fig. 2). A firm upper limit on the equivalent width of these lines is 1.7 mA which 
translates to a column density of 6x10^^ cm~^. At a temperature of 550 K, this corresponds 
to a total column density of C^^O of less than 1x10^^ cm~^. Allowing for the temperature 
uncertainty a 3a limit for ^^0/^''0 is >200. 



Discussion 



The basic characteristics of the HR 404 9 system are well un derstood. At the heart of HR 
4049 is a single-lined spectroscopic binary (IBakker et al.lll998l ). The visible early A/late B 
supergiant is a low mass, perhaps white-dwarf mass, post- AGE star. The unseen companion 
is an M dwarf or white dwarf of lower mass than the supe rgiant. The infra r ed pr ominent 
feature of the HR 4049 system is the circumbinary shell. lAntoniucci et al.l (120051 ) review 
the various geometries proposed for the circumbinary material. Con siderable evidence now 
points to a thick disk geometry. Detailed arguments are presented by lDominik et al.l (120031 ). 



In the following discussion we adopt the jPominik et al.l (|2003[ ) disk model (Figure 1) 
with the following key points. The disk is optically thick with a height-to-radius ratio ~l/3. 
The dust on the interior disk surface facing the star is approximately isothermal at 1150 K. 
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The temperature of the dust wall implies a distance between the star and dust of ~10 AU. 
The variability of HR 4049 and the hydrostatic scale height suggests that the inclination of 
the disk is ~60° (i.e. the plane of the disk is tipped 30° from the line of sight). The optical 
depth of the dust, the height of the disk, and the inclination result in only the far side of 
the disk being observable (Figure 1). 



5.1. Circumbinary Flow 



Previous observations of the CO first overtone are reported by iLambert et al.l (|l988h. 



Based on an excitation temperature of 300 ± 100 K and a non-stellar velocity iLambert et al. 



Jl988h conclude the CO is circumstellar. The much higher resolution and S/N data analyzed 
above refine the excitation temperature to 520 ± 20 K for the higher excitation lines and 
40±10 K for the low excitation lines. The velocities reported here and those reported by 
Lambert et al.l (119881 ) show no change over nearly 20 years, as expected for lines of circumbi- 
nary origin. Although of the current data is of higher precision, both data sets are consistent 
with an outflow velocity of ~1 km s^^. The column densities reveal that about four times 



more cool gas, ~2xl0^^ cm ^, is present than hot gas. 



-5x10^^ cm-2. 



The observations demonstrate that the gas is in rotational LTE and n ear or in vibrationa l 
LTE. For vibrational equilibrium the critical density is ylh ~ 10^° cm~'^ (iNajita et al.lll996l ). 
Taking this density and a CO column density of 2x10^^ cm^^, the thickness of the CO 
absorption line forming region is ~4xlO^^ cm. So the gas is restricted to a zone radially ~6 
R0 from the disk inner dust wall. The CO appears to depart slightly from vibrational LTE, 
so the density is likely slightly lower than the critical density. In any case, the thickness of 
the gas layer is certainly thin compared to the 2150 Rq spacing between the binary and dust 
wall. 

If the gas l ayer i s located radially just on the star side of the dust wall, adopting the 
Dominik et al.l (120031 ) geometry permits the total gas mass to be calculated. Taking the 
radius to be 10 AU and the height of the disk to be 1/3 the radius, the surface area of the 
cylindrical wall follows. The column density then gives the total number of CO molecules. 
Since the gas is oxygen rich, the number of CO molecule s is limited by the carbon abundance. 
Taking [C/H] for HR 4049 from IWaelkens et al.l (119961 ) and the solar carbon abundance of 



Grevesse et al.l ( 1l99ll ). the mass of the gas disk is 6 xlO^^ g m, i.e. ~ 0.1 M®. A total disk 



mass of a ^ 33 M0 was suggested by IPominik et al. J2003[ ). so the mass estimates are in 
accord with a thin gas zone at the edge of a more massive dust disk. 



The CO first overtone lines are symmetric and ~16 km s ^ across. In contrast, a line 
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width 20 times smaller, ~0.5 km s~^, is required to model the curve of growth. A 0.5 km 
s~^ width is consistent with thermal broadening. We suggest that the broadening to 16 km 
s~^ is due to a systemic flow of gas in the circumstellar shell. As noted above, our view of 
the dust disk continuum is limited to the side opposite from the star (Figure 1), so many 
kinds of axisymmetric flows could result in the observed line broadening. We consider the 
line shapes to constrain further our understanding of the flow. 

All the unblended 4.6 /im emission line proflles, including those for weaker lines, are 
double peaked. Since all the lines arc double peaked, wc discount self-absorption as the 
principal cause for this line shape. Double peaked lines suggest an origin in a rotating ring 
or disk. The observation of CO fundamental band absorption demands a P-Cygni type 
geometry, i.e., an emitting area extended relative to the continuum forming area. While 
emission is very dominant in the 4.6 //m HR 4049 spectrum, there is a hint of underlying 
and offset absorption in all the lines with the emission line profiles having a lower peak on the 
blue side than on the red side. The spectrum is dominated by saturated lines 10-20% above 
the continuum. However, some lines are stronger and we assume these stronger lines result 
as optical thick transitions cover larger areas. At 3 fim radiation from a 550 K blackbody is 
about half that at 4.6 fim and, indeed, the 3 OH Av=l lines are ~5% above continuum. 
Ultimately saturation combined with the physical extent of the line forming region demands 
a limited range of emission line strengths. The strongest lines in the 4.6 fxm CO spectrum 
of HR 4049 are ~40% above the dust continuum. 

The simple model of a rotating disk can be applied to the existing disk model (Figure 1) 
and tested by producing model proflles of the fines. The CO first overtone lines suggest that 
most of the CO occurs in a relatively thin layer. An absorption line was modeled by assuming 
a continuum source of 10 AU radius. The layer of absorbing gas was divided into zones of 
0.1 AU along the circumference. The fine RMS was assumed to be 0.5 km and the gas 
was assumed to be rotating in a Keplerian orbit. The resulting profile is a double peaked 
absorption line. This profile was then convolved with a Gaussian instrumental profile with 6 
km s~^ FWHM. The resulting synthetic profile, which is an excellent match to observations, 
is shown in Figure 13. 

Is this model also consistent with the 4.6 /xm emission line shapes? To investigate this 
question five assumptions were made: (1) The emissivity of the gas is constant over the 
entire region modeled. (2) The gas is in Keplerian orbits. (3) Line broadening is limited 
to the thermal b value, 0.5 km s~^ discussed above and the broadening from the Keplerian 
motion. (4) Absorption is insignificant. (5) The gas originates at 10 AU and extends to 
larger radii. Since the dust disk is opaque this extension is along the top of the disk (Figure 
1). An extension to larger radii was included since there is no requirement for a background 
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continuum source for the emission line spectrum. To fit the profiles with this model we found 
a maximum radius of 14 AU. As for the overtone model, the disk was divided into zones, 
the profile from each zone shifted and weighted by the viewing aspect, and then summed 
into velocity bins of 0.1 km s~^. The resultant double-peaked profile of the emission line is 
seen in Figure 14. Convolution with a Gaussian instrumental profile of FWHM = 6 km s~^ 
produced a good match to a typical emission line. 

While consistency between the modeled and observed profiles is satisfying, the funda- 
mental transitions clearly require much more refined modeling to address a number of details. 
For instance, there is a large differenc e between the 550 K C O and 1150 K dust temperatures. 



If, as is commonly assumed (see e.g. iGlassgold et al.ll2004l ). the gas and dust temperatures 



are in equilibrium within the disk then the 1150 K dust tempe rature applies only to a sur- 



mann 



19891) 



20M). It is 



face layer. Radiative cooling from CO fundamental emission flAvres fc Wiede: 
is largely disabled by the large optical depth of the CO lines (IGlassgold et al. 
plausible that the gas undergoes heating on exiting the disk. 

In an isothermal model optically thick CO self-absorption occurs for the fundamental 
transitions; the opacity in the low J 1-0 is ~400. The fundamental lines are seen in emission 
because the the 550 K temperature of the CO makes optically thick CO lines brighter at 4.6 
fim than the 1150 K continuum. At the resolution of the observations, 6 km s~^, narrow 
self-absorption lines of 0.5 km s~^ width are largely smeared out. Additionally, the gas 
is certainly not isothermal. If, for instance, the gas is heated as it leaves the disk, the 
temperature profile could increase toward the observer. Depending on the details of the 
spatial filling, optical depth, and temperature profile absorption is not a requirement. 

Two temperatures were measured in the CO first overtone, ~40 and ~550 K but no 
velocity differences were measured between the 40 K and 550 K regions, suggesting that 
these temperature regions are physically close together. Both the 40 K and 550 K CO are 
seen in absorption against the 1150 K continuum. H2O, on the other hand, is seen only in 
emission. Emission lines are not spatially limited to the 1150 K continuum forming region. 
If the absence of H2O absorption results from H2O existing only in the disk edge region and 
not in the disk mid-plane, the gas is differentiated vertically as well as horizontally relative 
to the plane of the disk. 

The measured decrease in the column density in between the v=0 and v=l levels implies 
that there is ample population to produce the observed optically thick 2-1 lines. However, the 
observation of optically thick 3-2 emission suggests that the v=3 level is populated above 
that expected from LTE. Overtone transitions higher than 3-1 are outside of the region 
observered. It would be of interest to search for the strongest lines in the 4-2 band. 
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The emission profile is shifted relative to tlie center-of-mass velocity by ~1.5 km s~^, 
suggesting an outflow. If the depression of the blue wing in the emission profiles is due 
to absorption in front of the dust disk, the absorption line profile is formed in a region 
with less outflow than the extended emission line forming region. Outflow was also noted 
for the flrst overtone CO lines. The outflow increases for the very lowest excitation lines, 
suggesting that the gas cools in the inner-disk region and is accelerated as it flows out. For 
the lowest excitation CO fundamental lines the cold outflow is seen i n absorption w i th the 
outflow velocity increasing (Figure 10) as excitation energy decreases. iDominik et al.l (120031 ) 
suggested that along the edges of the disk an outward flow results from ra diation pressure 
erosion. Alternatively, a disk pressure gradient can result in an outward flow (ITakeuchi fc Lin 



2OO2I ) without a need for small grains. Indeed, the outflow could be driven initially by either 



gas or dust since mora entum is transferred between the gas and dust through collisions 
jNetzer fc EhtzuilllQQsh . 



5.2. Comparison with Optical Spectra 



A detailed anal ysis of time series C I, Na I D (Di and D2), and H a spectra of HR 
4049 is presented by iBakker et al.l (I19981). The Na D lines contain a number of absorption 
components as well as weak emission. iBakker et al.l (119981 ) identify two Na D absorption 
components with the circumstellar environ ment of the binary system. These are labeled as 
'Ai' and 'A2' (see Table 2 and Figure 4 of IBakker et al.lll998l ). Ai has a velocity of ~ -5.0 
km s~^ (mean of Di and D2) relative to the systemic velocity. A2 has a velocity of -0.8 km 
s~^ again relative to the systemic velocity. A2 is stronger than Ai by about 50% and has a 
slightly greater FWHM. 

The continuum in the infrared is dominated by the dust continuum. However, in the 
optical the continuum is entirely from the stellar photosphere. Thus the optical absorption 
is formed along a pencil beam originating near the center of the circumbinary disk. The Ai 
velocity has similar velocity to the outflow seen in the lowest excitation 1-0 CO lines. This 
outflow is a cold wind perhaps leaving the system. The A2 outflow is close to the outflow 
velocity seen in the CO flrst overtone as well as the slightly higher excitation 1-0 lines (Figure 
10). This flow is an outward flow of warmer gas perhaps associated with circulation in the 
disk. 

Given that the star is the continuum source of the Na absorption and the rear inner 
walls of the disk are the continuum source for the CO absorption, perfect agreement is not 
expected in either line of sight velocity or in FWHM. The overtone CO has a much larger 
FWHM than the Na D, as expected given the larger range of velocities sampled by the CO 
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along the lines of slight to the CO continuum forming area (Fig. 13). 



Na I also has an emission component ('A3' in lBakker et al.lll998l ). This is perhaps due 
to fluorescent emission from the gas interior to the disk. The line proflle is disrupted by the 
Na D absorption components but the FWMH of the emission, ~21 km s~^, is comparable 
to that of the CO emission. 



The Ha line proflle is complex. iBakker et al.l (Il998l ) identifled two components, 'C^ax' 
and 'Rmm' which are stationary and presumably are associated with the circumbinary envi- 
ronment. Both are seen in absorption. Cmax has a large outward velocity, -21.3 km s~^. The 
velocity of Rmm is much less, -7.5 km s"^. The energetics of the Ha line are very different 
from those of the cold gas lines discussed in this paper. Ha also has an absorption feature 
'Bmm' which possibly varies in anti-phase with the primary. Detailed understanding of the 
excited gas sampled by Ha requires modeling beyond the current discussion. 



5.3. Properties of the Binary Members 



The conceptual picture of a thin gas layer co-rotating just in front of the dust wall 
suggests the observed velocities result from Keplerian rotation. A FWZl of 18 km s~^ 
implies a rotational velocity of 9 km s~^. Assuming Keplerian rotation and a 10 AU disk 
radius, the total b inary mass required is 0.9 Mq. This is in agreement with the total binary 
mass suggested by lBakker et al.l (119981 ). Bakker's mass was based on the mass function from 
the spectroscopic orbit and the assumption that the A supergiant had a typical white dwarf 



mass. 



The mass function from iBakker et al.l (119981 ). the total binary mass, and the orbital 
inclination allows a solution for the individual masses in the binary. We make the assumption, 
discussed in §5.6, that the binary orbit is co-planer with the disk. The deflnition of the mass 
function then yields the masses for the individual stars. The A supergiant has a very low 
mass of 0.58 Mq conflrmii ig the post-AGB stat e of this star. This mass, nearly equal to that 
of a typical white dwarf (IBergeron et al.lll992l ). implies that the mass-loss process for this 
star has terminated. The companion mass is 0.34 Mq. This mass does not resolve the status 
of the companion. While a mass of 0.34 Mq is low for a white dwarf and strongly suggests 
an M-dwarf, it is possible that the companion mass has been altered by evolution (§5.6). 



Lambert et al 



5.4. Winnowing 

(119881 ) report quantitative abundances for HR 4049 revealing an ex- 
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tremely met al-poor star with [Fe/H ] < -3 but near-solar C, N, and O: [C/H] = -0.2, [N/H]=0.0, 
[O/H]=-0.5. iLambert et al.l (Il988l ) argue that the ultra-low iron abundances found in a post- 
AGB star cai inot be primordia l since there are no known progenitor AGB stars with similar 
abundances. IVenn fc LambertI (Il990l ) and iBondl (Il99ll ) find similar abundance patterns to 
those in HR 4049 in the young main-sequence A Boo stars and gas in the interstellar medium 
(ISM). In all three cases the abundance pattern is deficient in refractory (high condensation 
temperature) elements but nearly solar in volatile (low condensation temperature) elements. 
This abundance pattern is explained in the ISM by the locking up of refractory elements in 
grains. 



Fi ye extremely iron-defic ient post- AGB stars are known in the HR 4049 class (Ivan Winckel et al 



19951 ). iLambert et al.l (119881 ) and iMathis &: LamersI (119921 ) have noted that all are A stars 
with no surface convection. A likely scenario is that the observed abundances result from 
peculiar abundances in little more than the observed photospheric layer. The very low re- 
fractory abundance in the HR 4049 stars results in a much lower opacity in the photospheric 
material th an from a normal comp osition making this region additionally s table against 



material tn an irom a normal comp osition maKing tnis region aaaitionaiiy s taple a) 
convection ( Mathis &: Lamers 19921 ). Assuming a 47 R© radius for HR 4049 (|Bakker 
19981 ) and referring to a 7500 K T^/f, log g = 1.0 model atmosphere (IKuruca 1 19791 ). the 



et al. 



photosphere of HR 4049 above optical depth unity contains a few percent of an Earth mass 
of volatile material. 



Mathis fc LamersI (119921 ) postulated that the HR 4049 abundance pattern results from 
the separation of mass-loss gas and d ust by differential forces on the gas and dust in a 
circumstellar shell. IWaters et al.l (119921 ) further suggested that a circumbinary disk played a 
critical role. Winnowing occurs as gas is accreted to the stellar surface while dust remains in 
the circumstellar shell or is ejected. The A Boo stars, which are also A stars without surface 
convection, have sir nilar surface abundance s due to winnowing of ga s from dust in a pre- 
main sequence disk (IVenn fc Lambertlll990l ). iMathis fc LamersI (119921 ) found the removal of 
refractory elements from a solar abundance gas to be a very inefficient winnowing process. 
They suggested that an efficient winnowing process is one that creates a gas with a low 
refractory abundance. 

Models of disks, created mainly to explore pre-main sequence evolution, provide a rich 
view of the basic disk physics. While the detailed physics of the winnowing are complex, 
these models, combined with the current observations, reveal the basics of the winnowing 
process. In optically thin circumstellar regions, the radiation pressure to gravity ratio for A 
stars dra gs grains with larger than 4 m i crons inward wh i le exp elling grains smaller than 4 
microns (ITakeuchi fc Artymowicg|200ll ). iTakeuchi fc LinI (120021 ) extend this to disk models, 
showing that large particles accumulate in the inner part of the disk. These models also 
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apply to opticall y thick disks where th e majority of the dust in the disk is not exposed to 
stellar radiation (ITakeuchi &: Lirul2003l ). In this case, interaction with the stellar radiation 
field at the inner disk edge drives flows in the disk with the dust-to-gas ratio increasing at 
the inner disk edge. 



Dominik et al.l (120031 ) conclude that the grain size distribution in the HR 4049 disk is 
currently indeterminate. However, they note that in the case where the inner disk consists of 
small grains, these grains will be driven outward by radiation pressure exposing a dust free 
region of gas. This gas layer will be driven inwards by either the gas pressure gradient or sub- 
Keplerian rotation. In pre-main sequence disks it has been shown that the gas interior to the 
dust suffers turbulent viscosity and accretes onto the central star. The visc ous time scale in 
pre-main sequence circumst ellar disks is typically estimated at ~ 10^ years (ITakeuchi fc Lin 
2OO3I : iHartmann et al.lEoQsh . 



The observations reported here of a sheet of gas at the inner disk edge support the 
model where separation of volatiles from grains occurs near the inner dust disk surface. The 
temperature of the gas released from the grains is far to low for evaporation of refractory 
elements to take place. The total gas mass interior to the HR 4049 dust disk is ~0.1 M^. 
The gas required in the stellar photosphere to alter the observed stellar abundances is one 
tenth this. A naive interpretation is that the surface material required to match the observed 
abundances could be accreted in ^ 10^ years. For post-AGB evolution this may be too long. 
An alternative is that the winnowing process currently observed is the termination of a very 
rapid clearing of the inner disk region that result in sudden accretion of the gas now present 
in the stellar photosphere. The observed CO, H2O, OH is near the disk. The flow to the 
stellar surface is presumed much more tenuous and is not observed. 



Takeuchi fc LinI (120021 ) found that higher than a few disk scale heights from the disk 



midplane the gas rotates faster than the particles due to an in ward pressure gradient . This 
drag causes particles to move outward in the radial direction. ITakeuchi fc LinI (120031 ) spec- 
ulate that in an optically thick disk, particles in the irradiated surface layer move outward, 
while beneath the surface layer, particles move inward. The outward flow seen in CO plus the 
driving entrained particles could rejoin the cool outer portions of the disk. In this case, the 
inward interior disk flow would move this material to the inner disk surface. The winnowing 
process could then be a distillation process resulting in a disk with increasingly refractory 
grains. 
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5.5. Isotopic Abundances 



The oxygen isotopic ratios reported by ICami fc Yamamural (l200ll ) set HR 4049 apart 
as having by a factor >10 the smallest ratios of ^^O/^^O and ^^0/^*0 known at that time. 
Our analysis of the optically thin CO first overtone transition does not support these results. 
There are no detectable C^^O first overtone lines giving a 3a limit of ^^0/^''0 > 200. On 
the other hand, the fundamental spectrum of CO in HR 4049 consists of optically thick 
emission lines. Four isotopic variants (^^C^^O, ^^C^^O, ^^C^'^0, and ^^C^^O) can be seen 
in the fundamental spectr um with lines of similar intensity. The CO2 bands measured by 
Cami fc Yamamural (120011 ) were observed at low resolution by ISO and are in the 13 - 17 fim 
region of the infrared. These bands appear in emission. We contend that the oxygen isotope 
ratios appear small in the se CO2 bands because, as f or the CO fundamental, the emission 
lines are highly saturated. ICami fc Yamamural (1200 ll ) warn that their isotopic ratios are in 
the optically thin limit. 



The carbon and oxygen isotopic ratios appear typical for an AGB star (lLambertlll988l ). 
While the oxygen isotopic ratio in the circumstellar environment of HR 4049 is not ab- 
normally low, there are stars that do have extreme oxygen isot ope values. Some hydroge n 
deficient carbon stars have ^^O/^^O considerably less than 1 (IClayton et al.l l2005l . 120071 ). 



Clayton et al.l (120071 ) suggest that the extreme overabundance of ^^O observed in these ob- 
jects is the result of He-burning in white dwarf mergers. Meteoritic samples have been found 
with small ^^0/^*^0 ratios. These could result from processes in the pre-solar nebula or 
pollution from a stellar source. While the rarity of stellar sources with small oxygen isotope 
ratios suggests a stellar source is unlikely, the ori gin of exotic oxyge n isotopic ratios detected 
in early solar system samples remains uncertain (lAleon et al.ll2005l ). 



5.6. Binary Evolution 



ISO data described by lDominik et al.l (120031 ) contain features from oxygen-rich molecules 
implying that the disk is oxygen-rich. The lack of mid- i nfrare d silicate features associated 
with oxygen-rich grains is attributed by iDominik et al.l (120031 ) to high optical depth in the 
dust disk. Our observations reveal a 2.3 - 4.6 /im spectrum resulting from a mix of gas 
phase molecules, CO, OH, and H2O, typical for an oxygen-rich environment. If, as seems 
probable, the gas consists of volatiles evaporated from grains then the disk environment is 
oxygen rich. 

In contrast, as reviewed in §1, carbon-rich circumstellar grains have been observed. An 
explanation is that these grains are exterior to the disk. Carbon-rich chemistry is the result 
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of evolution in the AGB phase where CNO material processed in the stellar interior is mixed 
to the surface. For AGB stars of mass ^ 2, the surface chemistry of the AGB star is converted 
to carbon-rich by the third dredge up. Rapid AGB mass loss then produces a carbon-rich 
circumstellar shell. In the case of HR 4049 the fossil carbon-rich shell of AGB mass loss is 
still observable although HR 4049 is now a post-AGB object. 



Why is the disk oxygen-rich? iBakker et al.l (119961 ) noted that the current binary separa- 
tion is less than the radius required for the AGB phase of the current post-AGB star. Hence, 
prior to the post-AGB stage the HR 4049 system underwent common envelope evolution. 
Prior to the common envelope phase the system passed through a pre-AGB contact binary 
phase with the more massive star transferring mass onto the less massive member. An AGB 
star does not contract due to mass loss, so the AGB star continued to expand enveloping 

the dwarf companio n. Common envelope systems rapidly eject mas s from both member s 

( iTaam fc Sandquistll2000l ). Most carbon stars have C/0 near unity (ILambert et al.lll986l ). 
Mixing or mass transfer during the common-envelope stage converted the carbon-rich enve- 
lope of the AGB star back to an oxygen-rich envelope. Mass lost during the common-envelope 
phase formed the current circumbinary disk. In such a sce nario a co- rotating circumbinary 
disk is formed surrounding the binary (IRasio fc Livid 1 19961 ). 



This process is apparently not unusual. iDe Ruyter et al.l (120061 ) find that circumbinary 
disks are a common feature of post-AGB stars. The current A-supergiant M-dwarf HR 4049 
binary is rapidly evolving to a white-dwarf M-dwarf binary system. White-dwarf M-dwarf 
binary systems with a co-binary di sk resulting from common-envelope evolution also appear 
to be common jHowell et al.lbood l 



6. Conclusions 



The 2 to 5 jjim spectrum of HR 4049 is formed in a circumbinary disk and wind. The 
optically thin 2.3 /im CO lines appear in absorption against the dust continuum, allowing 
the determination of the mass of gas. The gas forms a thin layer, of radial thickness ~ 6 
R0, lining the dust disk. This gas is composed of the volatiles separated in the disk from 
grains. The 4.6 /xm emission spectrum requires a region of line formation extended beyond 
the continuum forming region. The circumbinary gas is rotating with a Keplerian velocity 
of ~ 9 km s^^. Combined with the circumbinary disk radius and inclination derived from 
photometry, Keplerian rotation allows the determination of the masses of the individual 
binary stars. The very low mass of the A supergiant, 0.58 Mq, confirms the post-AGB 
nature of this object. Gas is also flowing out of the system, perhaps as a result of a disk 
pressure gradient, at ~1 km s~^. 
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Our observations show that the HR 4049 circumbinary disk has typical AGB abundances 
for the carbon and oxygen isotopes; ^^C /^^C = 6t^^ and ^^O/^^O >200. Exotic mechanisms, 



as proposed e.g. by lLugaro et aLl (120051 ). for pro duction of the °0/ "O are n ot required. The 



widely quoted value of ^^O/^^O ~ 8 reported by lCami &: Yamamural (120011) results from the 



naive interpretation that the infrared emission lines are optically thin. 



Cami fc Yamamura 



( I2OOII ) warn that their values were in the optically thin limit. The extreme saturation of lines 
in the 4.6 iim spectrum of HR 4049 results in nearly equal apparent strengths for isotopic 
variants of molecular species with abundances differing by factors of 10'^. 

The peculiar surface abundances of HR 4049 are likely the result of winnowing driven 
by the evaporation of volatiles in the disk and the viscous accretion of this gas onto the 
star. Detailed modeling of the process will be required to determine if there is adequate 
time for the current outgassing of the disk to fully explain the surfaces abundances of the A 
supergiant or if a sudden, post-common envelope clearing of the inner disk is required. The 
existence of the A Boo stars shows that the winnowing process applies to pre-main sequence 
as well as post-AGB systems. The wider significance of the winnowing process may well 
be in systems where the convective nature of the stellar photosphere cancels any impact on 
stellar abundances. However, circumstellar grains in these systems are undergoing processes 
separating volatile and refractory elements. This winnowing could have general application 
to the chemical evolution of grains in pre-main sequence disks. 

Carbon-rich circumstellar material implies that the post-AGB star was a carbon-rich 
star on the AGB. The current oxygen-rich circumstellar disk likely evolved from common 
envelope mixing. HR 4049 is one of five known post-AGB stars with similar photospheric 
abundances. Of the other four at least one, the red-rectangle nebula/binary HP 44179, has 



a sim ilar oxygen-rich circumbinary disk in a carbon-rich circumstellar shell (jWaters et al. 



19981 ). In future papers of this series we will explore the infrared spectra of the other 



members of the HR 4049 class of objects. 



This paper is based in part on observations obtained at the Gemini Observatory, which 
is operated by the Association of Universities for Research in Astronomy, Inc., under a coop- 
erative agreement with the NSF on behalf of the Gemini partnership: the National Science 
Foundation (United States), the Particle Physics and Astronomy Research Council (United 
Kingdom), the National Research Council (Canada), CONICYT (Chile), the Australian Re- 
search Council (Australia), CNPq (Brazil), and CONICRT (Argentina). The observations 
were obtained with the Phoenix infrared spectrograph, which was developed and is oper- 
ated by the National Optical Astronomy Observatory. The spectra were obtained as part 
of programs GS-2002A-DD-1, GS-2002B-DD-1, GS-2003A-DD-1, GS-2004A-DD-1, and GS- 
2005B-DD-1. We thank Drs. Claudia Winge and Bernadette Rodgers and the Gemini South 
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(JPL) funded by NASA through the Michclson Fellowship Program. JPL is managed for 
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Table 1. Log of Observations 



Date 


Wavelength 


Frequency 


S/N 






(cm~^) 




2002 Feb 13 


2.3120 


4324 


290 


zUUz reb io 


Z.OZOO 


4oUo 


ZZu 


2002 Feb 13 


2.3309 


4289 


280 


2002 Feb 13 


2.3421 


4268 


260 


2002 Feb 13 


2.3617 


4233 


400 


2002 Dec 11 


2.3405 


4272 


80 


2002 Dec 12 


2.3126 


4323 


90 


2002 Dec 12 


4.6434 


2153 


200 


2002 Dec 12 


4.6629 


2144 


190 


2002 Dec 14 


2.2980 


4350 


350 


2002 Dec 14 


2.9977 


3335 


130 


2002 Dec 14 


4.6219 


2163 


150 


2002 Dec 14 


4.6825 


2135 


170 


2003 Feb 16 


2.3416 


4269 


210 


2004 Apr 3 


4.8874 


2045 


220 


2005 Dec 10 


2.3634 


4230 


300 


2005 Dec 10 


2.3523 


4250 


350 
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Table 2. CO Av=l Line List 



Species 


Line 


i, 


1 


Line 


ic 


Line 


ic 


Line 


ic 


12^160 


1-0 


P3 


1.36: 


1-0 


P2 


1.25 


1-0 


PI 


1.27 


1-0 


RO 


1.28 




1-0 


Rl 


1 


.26 


1-0 


R2 


1.25 


1-0 


R3 


1.17 


1-0 


R4 


1.36 




1-0 


R5 


-1 

1, 


.38 






















2-1 


PIS 


-I 

1, 


.3d 


2-1 


R3 




2-1 


R4 


1.27 


2-1 


R6 






2-1 


R7 


1, 


.37 


2-1 


R8 


1.31 


2-1 


R8 


1.31 


2-1 


R9 


1.26 




2-1 


RIO 






2-1 


Rll 


1.27 


2-1 


R12 


1.27 


2-1 


R13 


1.22 




3-2 


PU 


i 

1, 


.13 


3-2 


RIO 




3-2 


Rll 


1.08 


3-2 


R12 






3-2 


R13 






3-2 


R14 


1.08 


3-2 


R15 




3-2 


R16 






3-2 


R18 






3-2 


R21 


1.09 
















5-4 


P5 


1 


.03 




















13^160 


1-0 


R9 






1-0 


RIO 


1.23 


1-0 


Rll 


1.20 


1-0 


R12 


1.23 




1-0 


R13 


1, 


.22 


1-0 


R14 




1-0 


R15 


1.27 


1-0 


R16 


1.20 




1-0 


R17 


1 


,20 


1-0 


R18 


1.21 


1-0 


R19 


1.17 


1-0 


R20 






2-1 


P7 


i 

1, 


.28 


2-1 


P6 


1.25 


2-1 


R17 


1.12 


2-1 


R18 






2-1 


R19 


1, 


.16 


2-1 


R20 




2-1 


R21 




2-1 


R22 


1.13 




2-1 


R23 






2-1 


R24 




2-1 


R25 


1.13 


2-1 


R26 






2-1 


R27 


1 


.09 


2-1 


R28 




2-1 


R29 




2-1 


R30 


1.08 


12^170 


1-0 


P17 


1 


.14 


1-0 


R3 




1-0 


R4 


1.11 


1-0 


R5 


1.15 




1-0 


R6 






1-0 


R8 




1-0 


R9 


1.12 


1-0 


RIO 






1-0 


Rll 


1, 


.08 


1-0 


R12 




1-0 


R13 












2-1 


R17 






2-1 


R19 
















12^180 


1-0 


P13 


1 


.26 


1-0 


P12 


1.24 


1-0 


RIO 


1.12 


1-0 


Rll 


1.14 




1-0 


R12 






1-0 


R13 




1-0 


R14 


1.12 


1-0 


R15 


1.12 




1-0 


R16 


1 


.18 


1-0 


R17 


1.13 


1-0 


R18 


1.11 


1-0 


R19 






1-0 


R21 


1 


.13 


1-0 


R22 


1.12 
















2-1 


R22 


1, 


.03 
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Note. — Central intensities at the maximum emission strength. Continuum 
= 1.0. 
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Table 3. H2O Line List 



Vibrational 


Rotational 


ic 


Rotational 


ic 


Rotational 


ic 


Transition 


Transition 




Transition 




Transition 




(lOO)-(OlO) 


[2,2,0]-[l,l,l] 


1.07 


[2,2,1]-[1,1,0] 


1.19: 


[3,l,3]-[2,0,2] 






[3,2,1]-[4,1,4] 


1.13 


[3,2,2]-[2,l,l] 




[4,0,4]-[3,l,3] 






[4,l,4]-[3,0,3] 


1.16 


[4,2,3]-[4,l,4] 


1.11 


[5,0,5]-[4,l,4] 


1.13 




[5,l,5]-[4,0,4] 


1.10 


[5,2,4]-[5,l,5] 




[5,4,2]-[5,3,3] 


1.03 




[5,5,0]-[5,4,l] 


1.06 


[5,5,l]-[5,4,2] 




[6,l,5]-[5,2,4] 


1.08 




[6,l,5]-[6,0,6] 




[6,2,5]-[6,l,6] 


1.07: 


[6,3,4]-[6,2,5] 






7,2,5-6,3,4 

L ' ' J L ' ' J 


1.05 


7,3,5-7,2,6 


1.07: 


8,2,6-8,1,7 


1.05 




[8,3,6]-[8,2,7] 




[9,2,7]-[9,l,8] 


1.08 


[9,3,6]-[8,4,5] 






[9,3,7]-[9,2,8] 




[10,4,7]-[10,3,8] 


1.03 






(OIO)-(OOO) 


[6,5,l]-[5,2,4] 




[8,8,1&0]-[7,7,0&1] 




[9,5,5]-[8,2,6] 






[9,8,1&2]-[8,7,2&1] 


1.04 


[10,3,7]-[10,0,10] 


1.02 


[10,3,8]-[9,0,9] 


1.04 




[10,4,7]-[9,1,8] 


1.06 


[11,3,8]-[10,2,9] 


1.08 


[11,7,4]-[10,6,5] 


1.04 




[11,7,5]-[10,6,4] 


1.03 


[12,4,8]-[11,3,9] 




[12,6,6]-[11,5,7] 


1.06 




[13,5,8]-[12,4,9] 


1.05 


[13,5,9]-[12,4,8] 


1.07 


[13,6,7]-[12,5,8] 






[14,6,9]-[13,5,8] 


1.06 


[15,6,10l-[14,5,9l 




[16,6,11]-[15,5,10] 




(OOI)-(OIO) 


[0,0,0]-[l,0,l] 




[1,1,0]-[1,1,1] 




[1,1,1]-[1,1,0] 


1.15 




[2,1,2]-[2,1,1] 




[2,2,0]-[2,2,l] 


1.18 


[2,2,0]-[3,0,3] 






ro 1 1 ro ni 
[2,2,iJ-[2,2,UJ 




Fq n q1 To nl 
[3,U,3J-[2,2,UJ 


1 on 








[3,2,2]-[3,2,l] 


1.23: 


[3,2,2]-[4,2,3] 


1.21 


[4,0,4]-[3,2,l] 






[4,2,2]-[4,2,3] 


1.06 


[5,0,5]-[4,2,2] 


1.02 


[5,1,4]-[4,3,1] 






[6,l,5]-[5,3,2] 


1.01: 


[6,3,3]-[6,3,4] 


1.02 


[7,l,6]-[6,3,3] 






[7,3,4]-[7,3,5] 


1.09 


[8,l,7]-[7,3,4] 




[8,2,6]-[7,4,3] 






[8,3,5]-[8,3,6] 


1.03 










(020)-(010) 


[5,5,0]-[4,2,3] 




[6,5,2]-[5,2,3] 


1.01 


[8,8,0&1]-[7,7,1&0] 


1.05 




[10,2,8]-[9,1,9] 




[10,3,8]-[9,0,9] 




[10,4,7]-[9,1,8] 


1.02 




[10,7,4]-[9,6,3] 


1.03 


[10,7,3]-[9,6,4] 




[11,3,8]-[10,2,9] 
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Fig. 1. — Cartoon of the 'wall' model for HR 4049 (top) taken from iDominik et al.l (120031 ). 
Below the Dominik model the spatially resolved observer's view of the system is shown in 
cross section. The disk is illuminated only from the inside. The cold disk blocks a large 
section of the inner 1150 K surface from view. The observer sees only that section of the 
1150 K disk inside the oval labeled "Cold Disk Vignetting Limit." 
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Fig. 2. — A selection of CO 2-0 R branch lines. The abscissa consists of ~8 A (1.5 cm^^ or 
~100 km s~^) increments of spectrum centered on each of the labeled lines. Top row shows 
i2q16q lingg^ middle row ^^CO hnes and bottom row C^^O lines. The columns aline the 
rotation quantum number J for the isotopic lines to approximately equal excitation (within 
J"±l). The spectral region containing ^^CO was not well covered, hence only a few lines 
are shown, however, ^^CO lines are clearly present in the spectrum. Only limits to the C^^O 
lines are detected. All CO first overtone lines have purely absorption profiles. 
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Fig. 3. — Overview of the 4.61 - 4.69 /xm spectrum of HR 4049 showing the forest of CO 
and H2O emission hnes. Gapped spectral regions indicate failure to restore the spectrum of 
HR 4049 due to optically thick telluric lines. 
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Fig. 5. — As per Figure 4 for the 4.64//m region. 
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Fig. 6. — As per Figure 4 for the 4.66/um region. 



-38- 




-cieo 1-0 jiia +(0 o i)[2 i 2]-[o i b)[z i i] 

13C0 2-1 R20 

-3-2 H13 



1-0 PI 

(1 0)[7 3 5]-(0 1 0)t7 2 6] 



O 

CV2 




170 1-0 RS 



13C0 1-0 HU 

(0 1)[0 0]-(C 1 0)tl 1] 
(0 1)[3 2 2]-(C 1 0)[3 2 1] 



00 

CO 

C\2 



-C180 1-0 R12 
13C0 2-1 R19 



(0 1 0)[6 5 1]-(C 0)[5 2 4]_ 




(0 1)[6 3 3]-(0 1 0)[6 3 4] 

-(0 1 0)[8 8 l,0]-(0 0)[7 7 0,1] _ 



-1-0 P2 

-<1 0)[4 4]-(0 1 0)[3 1 3] 



-2-1 R4 

(1 0)[10 4 7]-(0 1 0)[10 3 8] 

170 1-0 R4 
13C0 1-0 RIO 



13C0 2-1 RIB 
180 1-0 RU 



-3-2 Rll 

(1 0)[6 2 6]-(C 1 0)[8 1 7] 



_2_j-^(0 2 0)[10 2 8]-(0 1 0)[9 1 9] 

1-0 P3+(l 0)[6 3 4]-(0 1 0)[6 2 5] 

■(1 0)r3 1 3l-(0 1 0)[2 2] 

(10 0)[S 4 2]-(C 1 0)[S 3 3] 

13C0 1-0 R9 
170 1-0 R3 

13C0 2-1 R17 

180 1-0 RIO 

(1 0)[5 5 1]-(C 1 0)[5 4 2] 

(1 0)[5 5 0]-(0 1 0)[5 4 1] 
(0 2 0)[11 3 a]-(0 1 0)[10 2 9] 

(1 0)[4 2 3]-(0 1 0)[4 1 4] 
-3-2 RIO 



CO ^ 



CV2 



O 



CO o 

^ > 

02 CO 



CV2 
CO 



o 

CO 

cv 



C\3 



00 



CO 



Fig. 7. — As per Figure 4 for the 4.68/um region. 
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Fig. 8. — As per Figure 4 for the 4.89/um region. 
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Fig. 9. — An enlarged view of the spectrum shown in Figures 4 and 5 showing the regions 
surrounding the ^^C^^O R 2 hne (left) and the R 5 line (right). The R 2 line is unblended 
on the red wing while the R 5 line is unblended on the blue wing. Both lines have P-Cygni 
type profiles. Higher excitation CO lines as well as H2O lines shown in this Figure exhibit 
typical double peaked emission profiles. 
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Fig. 10. — Radial velocities of the absorption component of the CO 1-0 low excitation lines 
(J"=0 through 5) as a function of excitation energy of the lower level. A number of these 
lines are blended with other circumstellar lines. The bar to the right labeled ' high-J" ' is at 
the value of the mid-emission absorption for higher excitation lines. There is a clear trend 
for the low est excitation lines to have a larger outflow. The dashed line is the binary system 
7- velocity (IBakker et al.lll998l ). 
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Fig. 11. — The 3//m region spectrum of HR 4049 showing the OH 1-0 P2/5.5 and 1-0 P2e5.5 
hnes. 
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Fig. 12. — Boltzman plot for HR 4049 first overtone CO lines. Data are shown for the two 
isotopic species of CO that were detected in the first overtone spectra, ^^C^^O and ^^C^^O. 
The ^^C^^O excitation temperature is 40±10 K the low excitation lines and 530±20 K for 
the high excitation 2-0 lines. The four 3-1 ^^C^^O lines populate the 5000-6000 K region 
of the abscissa with an excitation temperature of 540±80K, suggesting a slight departure 
from vibrational LTE. However, a fit through all higher excitation lines remains within the 
uncertainties and gives an excitation temperature of 620±20 K. The ^^CO lines have an 
excitation temperature of 570±40 K. 
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Fig. 13. — Synthetic line profiles for the CO first overtone. The dot-dash line results from 
modeling a thin gas layer on the interior surface of the dust disk (see text). The dash line is 
the same model spectral line convolved to the instrumental resolution. The solid line is the 
observed profile of the CO R6 line. 
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Fig. 14. — A synthetic line profile for the CO fundamental lines compared to an observed 
profile. The dot-dash line is the synthetic line profile from an emitting zone near the 'cold 
disk vignetting limit' in Figure 1 (see text). The dash line shows this profile convolved to the 
instrumental resolution. The upper solid line is an observed CO line profile and the lower 
solid line is the difference between the model and observed line profile. 



